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ABSTRACT 

The continuous-time regularization scheme for defining phase-space path integrals is 
briefly reviewed as a method to define a quantization procedure that is completely covariant 
under all smooth canonical coordinate transformations. As an illustration of this method, 
a limited set of transformations is discussed that have an image in the set of the usual 
squeezed states. It is noteworthy that even this limited set of transformations offers new 
possibilities for stationary phase approximations to quantum mechanical propagators. 

1. INTRODUCTION 

For many years now it has been customary to define path integrals with the aid of 
coherent states [1]. Such formulations have been developed not only for the canonical 
coherent states suitable for the Weyl group (i.e., the Heisenberg algebra), but for coherent 
states based on other groups as well, notably the unitary and orthogonal groups with 
(in)definite signature, the affine group, etc. However, for the sake of convenience and to 
focus on the relation with standard squeezed states, attention in this paper will be confined 
to path integrals constructed with the aid of canonical coherent states. The construction 
of coherent state path integrals is generally carried out in one of two standard procedures 
[2]. To illustrate these two procedures let us first introduce a few standard definitions 
involving coherent states [2]: 

1 = J \p,q){p,q\dpdq/2TT, 

H{p,q ) = (j>,q\H\p,q ), 

H (p2,q-2-,pi,qi) = (p2,q2\H\pi,qi)/(P2,q2\pi,qi)i 
'H = Jh ( p,q ) \p,q}{p,q\dpdq/2Tr, 
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where \p,q), (p,q) € JR 2 , denotes one of a collection of coherent states defined by 

M) = e~ tqP e tp Q\0), [Q,P] = i , ( Q + iP) |0) - 0, 

all of which are normalized, (p, q\p, q) = 1. In addition, we have introduced two “symbols” 
associated with a fairly general operator 7Y, namely H and h as functions on phase space. 
In terms of these quantities, the propagator form time i' to time t" = t' + T, T > 0, is 
given by either of the two expressions 

n=N n=N 

]^[ (Pn+l,q n +l\pn,qn) e~ teH ^ n +^ qn + i; P n ^ JJ dp u dq n / 2 tt. 

72=0 7 ? = 1 


jf f\ 


\p\q') = 


N 


oo 


n=N 


/-/ 


n=N 


II (Pn+l,<ln+l\Pn,qn} J[ e ieh ^ Pn ^dp n dq n / 2tt, 

0 n=l 

where we have introduced the notation £ = T/(N+ 1), p" ,q” = p N+h q N+l , 
and p ,q = po,qQ. In a formal limit, in which the order of integration and the limit 
are interchanged and the integrand is evaluated for continuous and differential paths, the 
formal result emerges, respectively, that 


{p’\q"\e- mT \p\q') = M [Je^-^^VpVq, 

(p",q"\e- >nT \p',q') = M [ fe* S^-^M^VpVq, 

where, as conventional, we use a single standard integral sign here to represent a (formal) 
functional integration. Since, in the general case, H ( p,q ) ^ h(p. q ), we are seemingly 
led to a paradox, namely that two generally different expressions can be given for the same 
quantity. That these two expressions are different is just a dramatic reflection of the very 
formal nature of such equations in the first place; each is correct if interpreted in the 
manner indicated in the lattice regularized form given above. 


In recent years, a very different regularization and formulation of coherent state path 
integrals has been developed that is both rigorous in construction and does not exhibit the 
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paradox outlined above [3]. In this formulation, a continuous-time regularization scheme 
is found that takes the form 





where pf*, denotes a planar two-dimensional Wiener measure with diffusion constant 

v that is pinned so that p{t'),q(t') = p',q' ^ P (O > 9 (O = p" ^ * ^^ nce Ptf) 
and q (t) are (independent) Brownian motion paths, the integral Jp{t)dq(t) is properly 
understood as a well-defined stochastic integral [4]. In the present form the Ito and 
Stratonovich formulations yield the same result; however, under coordinate transformations 
the Stratonovich form is chosen and so it is convenient to adopt the Stratonovich form 
from the outset. It must be appreciated that the expression above involving the Wiener 
measure is rigorous and unambiguous; the marvel is that this genuine, i.e., continuous time, 
path integral formulation actually provides the correct propagator for the Hamiltonian 
7i provided one adopts the symbol h (p, q) to use as the classical Hamiltonian in the 
action even though it may, in general, contain a nonzero Ti. In a formal, but nevertheless 
suggestive language, one may also say that 





e * I\p<i-KPA)] dt e -to f[p 2+ v 2 } dt T>pVq, 


which shows the continuous-time regulatory nature of the indicated expression inasmuch 
as the v— dependent factor in the integrand formally goes to unity as v — ► oo. Although 
formal in nature, the last equation may be understood as a short hand expression for the 
former one when it is accepted that the various terms do not have independent meaning 
but only in combination with one another. Thus they may be recombined into the proper 
mathematical form at any time. (This is similar to how the “quotient” dy/dx should be 
understood for the derivative.) 


The expression for the propagator given above is not only well defined mathematically, 
but it also enjoys a covariance under generally time-dependent canonical coordinate 
transformations. Let two canonical coordinate systems be related according to the equation 

p dq — h (p, q) = p dq -\-dG ( p,q , t) — k (p, q,t ) 

that not only holds in the classical case but for Brownian motion paths as well thanks to the 
choice of the Stratonovich rule [4]. Under the same canonical coordinate transformation, 
the metric on flat space that supports the two-dimensional Brownian motion changes to 

dp 2 + dq 2 = do 2 {p,q, t ) = A (p,q,t) dpdp+B (p,q,t) dpdq -\-C {p,q,t) dqdq 
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and as a consequence the propagator takes on the form in the new coordinates, which for 
convenience we relabel p. q again, 

llm MtW- a ') I 

v — > oc J 11 

This formula expresses — for the first time and after over sixty years of the theory of 
quantum mechanics — a fully canonically coordinate covariant formulation of the process 
of quantization [5]. This expression clarifies the role of the Schroedinger quantization rule 
in which coordinates act as multiplication while momenta act as derivatives; this rule of 
quantization is valid in and only in Cartesian coordinates as often noted, but Cartesian 
coordinates in phase space rather than in /j-space and <y-space separately as commonly 
stated [6]. Of course, the arena for classical mechanics resides in a symplectic manifold and 
it does not employ a (Riemannian) metric in its formulation. On the other hand, quantum 
mechanics has a different and richer basis in which a metric structure appears. Indeed, 
it is not unreasonable from a classical viewpoint that a metric structure is appended to 
the classical phase-space manifold, not for purposes of defining the Hamiltonian equations 
of motion, but rather to keep track of just what physics a given system refers to. For 
example, an harmonic oscillator (centered at the origin) appears as an harmonic oscillator, 
e.g., with a Hamiltonian given by \ ( up 2 + 2 hpq + c.q 2 ) , a > 0 ; h > 0 ; ac. > b 2 , only in 
Cartesian coordinates in phase space. In non-Cartesian coordinates an harmonic oscillator 
assumes a different form from that indicated. Just what system actually corresponds to an 
harmonic oscillator (or free particle, or quartic anharmonic oscillator, etc.) is coded into 
the classical scheme by the implicit use of an auxiliary flat metric on the two-dimensional 
phase space, and its expression in Cartesian coordinates. This same flat metric space 
actually enters the formulation of the quantization procedure as described in the present 
article through its use as a carrier for the Brownian motion. Once it is decided which sets 
of canonical coordinates are the Cartesian ones, so that the expression for a system which 
represents (say) an harmonic oscillator is unambiguous, then the quantization procedure 
itself is unambiguous in the approach advocated here. After the well-defined path integral 
is set up, then one is free to make a variety of coordinate changes within that integral, 
among which possibly time dependent canonical transformations are to be distinguished. 
Indeed, one can go so far as to introduce a Hamilton-Jacobi transformation so that the new 
Hamiltonian vanishes. This puts all the dynamics into the curvilinear coordinate system 
that is used to track the two-dimensional planar Brownian motion. As a consequence, the 
overall level of difficulties is conserved, as one would expect to be the case. 

It is hard to illustrate this program in its full potential, but it can be shown in a sort 
of small scale fashion. Indeed, it is squeezed states that can be used to provide a limited 
illustration of this overall program, and it is this “miniature” illustration to which we now 


250 



turn our attention. A convenient place to start the investigation is with the kinematics 
rather than the dynamics, and this in turn can be done simply by looking at the propagator 
for vanishing Hamiltonian. 


2. CHANGE OF VARIABLES IN THE PATH INTEGRAL: 

CONSTANT w TRANSFORMATIONS 

KINEMATICS 

If the Hamiltonian vanishes, or in the limit that T — ► 0, the “propagator reduces to the 
reproducing kernel , an integral kernel representing a projection operator onto the relevant 
subspace of all square integrable functions on phase space as given by 

l|pW;l> = 


which may be explicitly evaluated as 


IbW; 1) = jWrtu’-o-ilv-rtw-'rr]. 


In these expressions we have added a “1” to the label to emphasize that the coherent states 
are those based on an harmonic oscillator ground state with a unit angular frequency, 
w = 1. In particular, for a general value of cu, the configuration space representation of 

the coherent states reads 


( t J 
- 


-i u)(x-q) 2 +ip(x-q) 


and it follows that the overlap of two such states for the same value of w is given by 


(p",q"-,u\p',q';u) = J(p",q";u\x){x\p' 1 q'-,u)dx 

= exp {i (p" + p') (?" -?')-? (?" - P') 2 + («" 

It is clear therefore that the coherent state overlap obeys the identity 




which relates dilation of the angular frequency to a corresponding dilation of the coherent 
state labels, i.e., an expansion of one phase space coordinate and a contraction of the other. 


This relation may be codified another way as well. Let 

D = \ {PQ + QP) 
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denote the self-adjoint dilation operator with commutation properties [Q,D] = 

iQ, [P,D] = -iP, then it follows that coherent states for different angular frequen- 
cies are connected by 

| p,q\u) = e l 2 ln(a,)z? |cu - 2p,u;2g;l). 

Thus the unitary transformation generated by D is nothing other than the squeeze operator 
relating coherent states and squeezed states, or relating two sets of squeezed states with 
different squeezing values. In forming the overlap illustrated above, the squeezing operator 
drops out leading to the indicated relation. 

A path integral expression for the coherent state overlap at angular frequency lj can 
be readily obtained just by a coordinate change of the path integral appropriate for a unit 
angular frequency. In particular, if one makes the change of integration variables given by 

p(t) -*■ q (t) -*■ u*q (t) , 

then it immediately follows that 

(p",q";u\p',q'-,u) = „ J l^'F+^WVpVq 

showing quite clearly the connection of the relative scale factor in the two-dimensional 
Brownian motion and the parametric dependence in the coherent state representation. All 
this has assumed that u has been constant throughout; next we take up the case of a time 
variable u>. 


3. CHANGE OF VARIABLES IN THE PATH INTEGRAL: 
NONCONSTANT TRANSFORMATIONS 


The overlap of two coherent states for two different values of w is given by 


ip" ,q"-,u”\p' ,q'',u') = J {p" ,q"\J'\x)(x\p' ,q'\J)dx 


V 2 


\J 'fil? 


exp 


. {p"J + p'u") {q" - q') (p” -p ') 2 


u" + w' 


2{u" + cj f ) 


J'uj' (q" - q'f 
2 (c o" + J) 
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This expression also exhibits an alternative form given by 




where we have introduced a smooth but otherwise arbitrary function u; ( t ) t < t < 
t", which interpolates between u" = oj{t") and a/ = u>(t'). Of course, this expression 
holds as well even in the special case that u> ,f = u/ in which case u (t) goes smoothly 
between equal initial and final values, but is otherwise arbitrary. 


In the latter form the overlap of two coherent states for differing angular frequencies 
has been expressed in terms of the matrix element of a kind of propagator between coherent 
states of the same angular frequency. But the latter form admits a path integral expression. 
In particular, it follows that 


(p",q";l\e- i '% Ddt \p , ,q , ;l) 

= hm M 
V ► oc J 

Now, much as was the case earlier when u> was constant, we next make a time-dependent 
change of variables of the form 


p(t)-> w(<) 
q{t) -* u(t)iq(t), 


where u" = uj{t ") , <J = u) (t f ) . In making a time-dependent substitution of variables, 
additional terms will arise in the path integral integrand on the right hand side. In particular, 
the term 

S\pq- £ p<i] dt - 

the formal flat measure remains unchanged, 

VpVq — » VpVq, 

while the all important formal weighting factor 

e ~~h /[^ 2 +9 2 ]^ _> e ~b /[ ta, (0 _, p J + u, (09 2 ]^ t 
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Other terms might be contemplated in the exponent of the final expression such as those 
involving time derivatives of uj (t)\ however all of these will be negligible in the limit that 
v — > oc since they are not as singular as the indicated terms. While we prefer this heuristic 
characterization of the transformed Wiener process one should bear in mind that only the 
coordinate description of the planar two-dimensional Brownian motion is being changed 
and the process itself is in no way effected. We are encoding this change of coordinates 
by means of the change of coordinates of the metric on the plane. (A rigorous analysis of 
such transformed Wiener processes is in progress by the present authors.) Thus it follows 
after such a substitution of variables that 


M J e c ~ ± -f H ' ) ' P 2 +"(f VpVq. 


Consequently, the introduction of a smooth, time-dependent angular frequency that 
interpolates between the initial and final values in the Wiener measure provides just the 
right ingredient to yield the overlap between two different coherent states based on two 
different angular frequencies. This expression yields a simple but nontheless bona fide 
example of how the classical Hamiltonian — here just u (t)p (/) q ( t ) /2a? (t) — may be 
eliminated in favor of a change of coordinates with which to describe the two-dimensional 
Brownian motion on the phase space plane. Such an elimination additionally involves a 
change of coordinates at the endpoints, as illustrated in the central equation, but in the case 
of squeezed states, there is an alternative interpretation involving coherent states based on 
differing angular frequencies as embodied in the first part of the equation. We now turn 
our attention to the inclusion of dynamics in this example through the presence of a rather 
general nonvanishing Hamiltonian. 

4. CHANGE OF VARIABLES IN THE PATH INTEGRAL: NONCONSTANT 
(j TRANSFORMATIONS AND A GENERAL HAMILTONIAN 

INTRODUCTION OF GENERAL DYNAMICS 

Based on the earlier discussion it is quite straightforward to include a rather general 
Hamiltonian h ( p , q) into the problem. In particular, based on the initial discussion, let 
us consider 




v — ► oc 


M J e* Hf+fWVpVq 
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which after a change of integration variables becomes 


lim 


v 


oo 


where 


M 'p'+“'(>)?] d, T>pVq 

= (u/'-Vv'V; i|r e -v[«+* D l ,l, |u'-V,w'V; i) 

= (o/'-Vv'V; i| e -t'"«"B re -.7w'(*M< e ita<-'fl| t /-iy )W 'i 9 ' ; i) 

= <p", «"|Te-^ W ' (,)d, |p', <A 

fc (p (<) , 9 (<) , <) = ft (« (t)~‘ p (t) , w (t)iq (tj) , 


which contains an explicit time dependence from the angular frequency as well as an 
implicit dependence just from the time dependence of p and q themselves, and in addition 

where 

H'{t) = e*» lu ln ~^ D . 


The basic significance of the preceding equations can be summarized as follows: 
(p",</’^"\Te-‘l n 'W\p' t q';uj') 

= <</'- VV'V; i|Te- i ^ w+ * £ ’]'“|w'- VW; i> 


lim 


v — * oo 


M J e* /[P9- fc (P»9)“ ttP a +?} dt VpVq 


p(t)-»u/(t) ^p(t) 

1 


lim 


V — ► oo 


l P^ u ^] dt T>pVq. 


This is the most general form we are able to offer using squeezed states, and it shows, in the 
first of the equalities, how part of the Hamiltonian can be absorbed quantum mechanically 
by a change of the fiducial vectors — indeed just like going to the interaction picture in 
ordinary quantum mechanics, which is then responsible for the introduction of the time- 
dependent “interaction” picture Hamiltonian H' (<). The second pair of equalities just uses 
the original form of the path integral as modified by a change of variables that effects 
the end point conditions as well. The final equality just accounts for that very change of 
variables as requested in the line above. 
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5. EQUAL END POINT ANGULAR FREQUENCIES 

Let us return to the path integrals discussed at the beginning of this article, namely 
to those for which the initial and final angular frequencies are the same. For the sake 
of convenience, let us choose that value to be unity, i.e., J' — J = 1, and return to 
the original notation for the coherent states with unit a ;, namely that j p,q) = \p,q-, 1). 
However, this time we will retain the option of using a time-dependent angular frequency 
to (t) to interpolate smoothly between the original and final values of unity. In this case 
the formulas developed above simplify to become 


lim 

v — ► oo 


P( 0“^(0 2p ( t ) 


lim 


M J e i }[pq-h(p,q)--%.pq]dt e -±; J\p*+q*]dt J)pVq^ 

9(0-w(A(0 

M J e* f\Pi- k (PAJ)]<it e —z; J[u(t)~ l p : ‘+u(t)q*]dt Dy'Dq 


u — ► oo 

With a slight generalization, we can now turn this equation around to read 

jf J 


(P,q"\ e ~ ini bW> 


lim 

v — ► oo 
lim 


v — ► oo 


Mje i}{pq-h(p,q)}dt e -± j\p 3 +q*]dt X>pVq. 

M J S\p^q Q )di p) p Vq 


pUJ-MO ^p(t) 
1 

$( 0-^(0 2 q ( t ) . 


lim 


is — > oo 


M / l p 3 +u{t)q a ]dt x>ppq. 


In this expression the third line holds because it simply corresponds to a change of 
integration variables that does not have any effect on the values of the boundary labels 
since us" = J = 1. The last line represents, just as before, the consequences of that very 
change of variables. 
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Now observe that on the left side of this equation there is no reference to the function 
w (<), t' < t < t", u(t") = u(t') = 1, while on the right side of the equation, in 
the last part of the equation in particular, the function to ( t ) enters in a prominent way. 
This becomes especially significant when an approximate evaluation of the path integral 
is admitted, such as that which arises from a stationary phase approximation. Stationary 
phase approximations for coherent state path integrals with Wiener measure regularization 
of the kind considered here have been worked out previously [7] and we do not repeat 
that discussion here. The point we wish to emphasize, however, is that the choice of the 
angular frequency u (t) mil enter most probably in the form of the approximate solution, 
and naturally some expressions will be better approximations to the real answer than other 
expressions will. Just which will be the best approximation is, of course, not too easy to 
establish. Perhaps one scheme is to ask that the result be stationary with respect to small 
changes of the functional form of co ( t ). In practice one might want to let u> ( t ) depend on 
just a few discrete parameters and to seek stationary variations with respect to just these 
few parameters. This certainly seems easier to do than to ask for an extremal variation 
with respect to the entire function u (t). 


One can actually see a miniature working of this general kind of procedure in comparing 
the usual and the Maslov stationary phase approximations to the sharp position propagator; 
see, e.g., [7]. As given earlier, the configuration-space form of the coherent states given by 


(x\p,q;u) = 


)* 


e -%u>(x-q) 2 +ip(x-q) 


makes clear that (x\p, q-,u) = 6(x-q) which converts the coher- 

ent state representation to the sharp position or configuration representation, while 
lim / l w / T \ n a . u \ _ 1 e ip{x-q) which converts it to the dual or momentum rep- 

resentation (up to an unimportant phase factor). These features can also be seen in the 


relation 


(jp" ,q"',oj"\p = J {p" ,q"',u"\x)(x\p' ,q ;u )clx 


V~2 



exp 


LU | / UJ 

Z7 i " V u" 


* + J 


(p"-p') 2 

2 K + u>) 


(/ - q'f 
2 (c o" + w') 


that gives the overlap of two coherent states based on differing angular frequencies. 
Consider the limiting situation in which both uj" — * oo, <J — > 0. In that case it follows that 


lim 

bj" — MX) 


167 V 2 CJ f J 


1 

4 


(p", 


1 Jp'W'-i') 
s/T* 
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Thus it should be no surprise that the two standard stationary-phase type approximations are 
actually contained in the coherent state approach in the form of suitable limits. In a manner 
of speaking, the usual configuration space approach just involves choosing a constant and 
very large value of the angular frequency parameter ui (taken to infinity at the end of the 
calculation) and making a stationary phase approximation to the resulting path integral. 
On the other hand, the Maslov approach takes the propagator from a sharp configuration 
initially to a sharp momentum finally, approximates that path integral by a stationary phase 
approximation, and then returns the end point to configuration space by a Fourier transform. 
This approach can be approximated in our method by taking an angular frequency history 
u ( t ) that is initially huge (tending toward infinity) and finally very small (tending toward 
zero), approximating that path integral by a stationary phase approximation, and finally 
making a change from coherent states based on a very small angular frequency to one based 
on a huge angular frequency just by the kinematical factor given above. The coherent 
state approximation developed in particular in reference [7] proceeds in yet another way, 
namely, starting with a sharp configuration initially, propagating to a coherent state with 
a finite nonzero value of the angular frequency, i.e., ui = 0(1), approximating that path 
integral by a stationary phase approximation, and then passing from the final coherent state 
representation to a sharp configuration one. This approach can also be approximated in our 
scheme by having an u> ( t ) that initially is huge, and finally is finite and nonzero [ 0(1) J, 
approximating, in turn, that path integral by a stationary phase approximation, and then 
passing back to a coherent state based on a huge angular frequency at the final point. 


6. CONCLUSIONS 

In this article we have attempted to show the reader what the authors believe is 
the “latest” in path integral construction — the state of the art — and illustrate how 
variable changes can be rigorously carried out within the path integral formulation itself. 
Squeezed coherent states have been used as convenient bases throughout in the illustration 
of the general program by a “miniature” subprogram involving a fairly limited change of 
integration variables. The resultant formalism is able to express a path integral in terms 
of an essentially arbitrary function, the time varying angular frequency, u (t), which lends 
itself to various selections in case an approximation scheme is invoked. By illustrating 
that the usual and the Maslov approaches are but two small examples of how such an 
optimization can be used, it becomes clear that there are hidden in these formulas a whole 
host of differing approximation schemes some of which, in certain applications at least, 
may well be better than the schemes currently in use. It is left to the future to see just 
how to exploit the vast number of possibilities that have been opened up here. 
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